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1. Introduction 



, ^Department of Physics and Astronomy, SUNY, Stony Brook, New York 11794, USA. 

^ , The results of a number of constituent quark models in matter may be understood in 

\^ I the mean-field approximation by using a simple four-fermi model in 0-1-1 dimensions. 

o 

00 ' The fate of the spontaneous breaking of chiral symmetry in matter is still stirring a 
^ ! broad interest in light of present and future heavy-ion collision experiments. Most of our 
' understanding of the issue from first principles is limited to finite temperature, where 
I I lattice simulations have been carried out [|1| . At finite chemical potential, present Monte- 
Carlo algorithms are upset by the complex character of the measure . 
^ ■ A number of past and present analyses of the subject relies on QCD inspired models, 
such as the instanton model 0,^, or variants of the NJL model 0]. A common feature 
^ [ to all these models, is the the emergence of chiral constituent quarks in the vacuum, 
^ I following the spontaneous breaking of chiral symmetry. These models do not confine. 

The purpose of this talk is to show that the underlying mechanisms at work in most of 
these models can be captured using a four-fermi model in 0-1-1 dimensions 0. In section 
1, we review the model in the presence of chiral scalar and vector interactions. In section 
2 we analyze its schematic thermodynamical content. In section 3, we construct the quark 
condensate and discuss its relation to the scalar quark content of the constituent quark at 
low densities. In section 4, we construct the resolvent and the distribution of eigenvalues 
for the present model in the quenched approximation. In section 5, we derive the spectral 
sum rules for the present model in the quenched and unquenched approximation. Our 
conclusions are given in section 6. 



2. Model 

Consider a quark field ipaj where a = 1, 2, iV are 'color' indices, and f = 1,2, Nf 
are flavor indices. For simplicity Nf = 1 unless specified otherwise. By analogy with NJL 
models, the Lagrangian density is chosen to be 



£i = tjj^i'y 484 + im + ifi'y 4)11) (1) 
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or equivalent ly 



+ L^l{^^ + {^i-u^))R (2) 
+ RH{P + m)R + LH{P^ + m)L . 



The model will be discussed on a circle of radius [3 = l/T in one dimension unless 
specified otherwise. The gamma matrices are chosen such that 74 = offdiag(l, 1) and 75 = 
diag(l, —1), in the chiral basis with ip = {R,L). The auxiliary fields are: P = —2ig'^L'^ L, 
P^ = -2ig^R^R and u;4 = -2c/2^t^^^^. 

For m = (chiral limit), the model exhibits massive (constituent) quark excitations at 
low densities, and massless (free) quark excitations at high densities in the limit N ^ 00 
(mean-field). In the same limit, qq, qqq, qq, ... phases and/or exitations are down by 
Although the model lacks confinement, it bears much in common with more realistic chiral 
constituent quark models such as the instanton or the NJL model 0]. For an analysis of 
the thermodynamics in a confining model we refer to 0. 

3. Thermodynamics 

On a circle with boundary condition tplr + (3) = —iIj{t), the operator id^ is invertible 
with a discrete spectrum ujn = (2^ + 1)ttT. In the mean-field approximation or large A^, 
the pressure per particle associated to (0-0) for A^j = 1 is 

-p = uj -Tlog(l-n)n-n) ~^PP^ + auJ4Ujl (3) 
n 

with the usual occupancies n,n = 1/(1 + exp((ct; =1= (/i — Ui))/T)). Here u = \P + m\, 
S = Yi/2g'^ and a = V^/g^, following the rescaling ip — * v^V- The thermodynamical 
limit will be carried with A^ ^ 00 at fixed n = N/V3. (^ describes a quark with two 
energy levels ±00 at finite temperature and chemical potential. The first term is the 'zero 
point' motion, and the last two terms are exchange contributions ||^. The gap equations 
are 

2SP = 1 — n — n , p = n{n — n) = 2an uj^ . (4) 

(^ admit several solutions of which the one with maximum pressure will be selected. 
Generically there are two phases: a broken phase with constituent quarks and a symmetric 
phase with free quarks. For a < E there is a /i region with no real solutions to @). In the 
symmetric phase the pressure is p= n (/x— l/4a), while in the broken phase it is p = n/4S, 
for m = T = 0. The phase change occurs at pc = 1/4S + l/4a. 

For a = 00 (vanishing vector coupling) the model is that used originally by us and 
others In the broken phase P* = 1/2S for p = with the pressure p = n(m + l/4E). 
In the symmetric phase P* = for p = n and p = np. The phase changes are mean-field 
driven and can be analyzed at the critical points using universality arguments 0. For 
any m there is a first order transition at p^, = (m-|-l/4S) 
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Figure 1. Isotherms at m = (left) and m = 1 (right). 



The occurrence of a first order transition for small values of the temperature for light 
current masses is not generic to constituent quark models. Indeed, in NJL models the 
first order transition in the chiral limit is usually turned to a cross-over transition for light 
quark masses ^ . The cross-over is robust against parameter changes in the presence of a 
vector interaction. The first order transition observed in these models at T = is much 
like the one studied in Walecka-type models in the form of a liquid-gas transition at low 
nucleonic matter density [Q. In the latter, the nucleon mass is the constituent mass and 
m = 0. 

In Fig. |1| we show the isotherms in temperature steps of AT = 1/2 for m = (left) 
and m = 1 (right), with E = 1. At low temperature the transition is first order with 
a jump in the density (Maxwell construction) turning into a second order transition at 
high temperature. This is in agreement with our original arguments 0,0, and more 
recent investigations |]TU[. With increasing mass the 'tricritical' point shifts down to lower 



temperatures, and the second order transition is turned to a cross-over. This confirms 
what is usually observed in NJL type models: a first order transition in the massless case 
and a cross-over in the presence of light quarks p[ . 



4. Quark Condensate 

The behavior of the quark condensate \{qq)\ = 2nEmQ with mq = and m = is 
shown in Fig. ^ as a function of temperature and density in the chiral limit. The solid line 
is the phase boundary crossing the tricritical point (big dot) beyond wich a mixed phase 
is developing. The hole in the middle reflects on the mixed phase. The size of the hole 
shrinks with increasing mass m, making it closer in shape to the result obtained in NJL 
type models in 4 dimensions Above the tricritical point T^c. = f^q/'^ and = n/v^ 



Figure 2. \{qq) \ in the Tp plane. See text. 



or /i3c = mg log(2 + A/3)/3, the transition is first order and disappears at /i* = mQ/2. At 
/i = 0, the transition is second order (mean-field) and sets in at X,, = mQ/2. Amusingly, 
for a constituent mass mg ~ 300 MeV, p^ = 150 MeV, and p^c ~ T^c ~ 100 MeV. 
In the absence of matter, n = ~ (200 MeV)^. Hence pzd Po ~ 10/3 where po = 0.17 

fm~^ is nuclear matter density. Most of these numbers are similar to the ones observed 
in NJL models [^. 

The slope of the chiral condensate for small densities carries information on the 'pion- 
nucleon' a term. Here the role of the nucleons is played by constituent quarks. For 
small temperature the slope is —1 at the origin and refiects on the scalar charge of the 
constituent quark. This is to be compared to —2.5 from the 2-fiavour pion-nucleon sigma 
term with do ~ 40 MeV and a current mass m ~ 8 MeV. 



5. Quark Spectrum 

For a = oo, the pseudoscalar four- Fermi interaction in (|ip causes the quarks to interact 
as if they were moving in a random Gaussian potential provided by the new auxiliary 
fields Aab{T) ~ i'a'^liT) which is an X complex valued function of Tm|. For fixed 
A, the quark spectrum follows from 

-i4{idi + A + ip)^^ = \k[A]^^ (5) 
with anti-periodic boundary conditions and 

2^ = A(l + 75) + At(l-75). (6) 

which is symmetric, complex and block-off-diagonal. (^ in the present model is the 
analogue of the QCD eigenvalue equation in external gauge field. The Gaussian averaging 
with moments 

{AabA^cd}^ = ^^SacSbd (7) 

restores the four-fermi interaction. The spectrum (in the n = sector for A symmetric) 
can be readily probed through the distribution of eigenvalues 

giz,z) = j^j:{5iz-X,[A]))^ (8) 

k 
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Figure 3. Spectral sum rules. See text. 



which is complex valued for finite /i. Setting T = in (^) the model dimensionally reduces 
to a 0+0 dimension one which is a matrix model. 



The quenched distribution (g) can be easily constructed |TT],|T2[. For small values of 
H the density of eigenvalues concentrates near the imaginary axis, while for large values 
it does not. The unquenched distribution of eigenvalues is harder to generate (complex 
weight induced by quark feedback). However, we expect qualitatively a similar behavior 
but with a finite accumulation of eigenvalues on the real axis. The chiral transition from 
finite density to zero set in at intermediate values of fi for which the density of eigenvalues 
along the real axis would vanish in the thermodynamical limit N ^ oo. To characterize 
the onset of the transition quantitatively we suggest to use the microscopic sum rules as 
we now show. 



6. Spectral Sum Rules 

For fixed Nf, the spectral sum rule for our case reads 

where the averaging is over the matrices A with the unquenched measure in the 0+0 
model (no Matsubara modes). The sum rule (^ allows for a determination of the chiral 
condensate (low density). Near // ~ yU^, the finite size corrections are important, upsetting 
the present construction. In Fig. ^ the chiral condensate extracted from (^ in the 0+0 
model (averaged over 5 x 10^ matrices) is shown for one (upper left) and two (lower left) 
fiavours. In the upper right corner we show the chiral condensate extracted from the 
next sum rule. The dashed lines indicate the analytical result of the condensate. The 
lower right corner shows the sum rule (|^) for the quenched [Np = 0) case and finite N 
(for oo (H) diverges). Due to the elongation of the support of eigenvalues for large 

chemical potential, the sum rule changes sign. The numerical analysis indicates that this 
change occurs around the critical chemical potential (/x* ~ 0.5277 |]lll). This suggests 




(9) 
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that the quenched physics may still 'remember' the unquenched one at the critical point, 
through the way the complex eigenvalues in (|^) get redistributed near zero. 

7. Conclusions 

We have constructed a chiral four fermi model with the thermodynamical structure of 
a two-level quark system. This model undergoes first and second order-type transitions, 
and can be used to illustrate how the quark spectrum may be probed in matter. We 
used this model to show that the microscopic spectral sum rules are sensitive to tem- 
perature and chemical potential changes including a phase transition in the unquenched 
approximation. The interplay between the thermodynamical limit, the chiral limit, the 
quenched-unquenched approximations and the precision of numerical algorithms are read- 
ily addressed in this model. We believe the present analysis to be useful for current lattice 
QCD simulations with matter. 



8. Acknowledgments 

GP and IZ thank Gerry Brown and Madappa Prakash for discussions. This work was 
supported in part by the US DOE grant DE-FG-88ER40388, by the Polish Government 
Project (KBN) grants 2P03B04412, 2P03B00814 and 2P03B08614 and by the Hungarian 
grants FKFP-0126/1997 and OTKA-T022931. 



REFERENCES 



1. C. Detar, [hep-ph/ 9504325| ; and references therein. 



2. J.B. Kogut, M.P. Lombardo and D.K. Sinclair, Phys. Rev. D51 (1995) 1282. 

3. D. Diakonov, |hep-ph/9602375| ; T. Schafer and E.V. Shuryak, Rev.Mod.Phys. 70 



(1998) 323. 

4. M. Nowak, M. Rho and I. Zahed, 'Chiral Nuclear Dynamics', World Scientific 1996; 
and references therein. 

5. for a review see U. Vogel and W. Weise, Prog. Part, and Nucl. Phys. 27 (1991) 195; 
S.P. Klevansky, Rev. Mod. Phys. 64 (1992) 649; M.K. Volkov, Phys. Part. Nucl. 24 
(1993) 35; T. Hatsuda and T. Kunihiro, Phys. Rep. 247 (1994) 221. 

6. R.A. Janik, M.A. Nowak and I. Zahed, Phys. Lett. B392 (1997) 155; R.A. Janik, M. 
A. Nowak, G. Papp and I. Zahed, Acta Phys. Pol. B27 (1996) 3271 and Acta Physica 
Polonica B28 (1997) 2949. 

7. T. Hansson and I. Zahed, Phys. Lett. B309 (1993) 385, and references therein. 

8. T. Wettig, T. Guhr, A. Schafer and H. Weidenmiiller, |hep-ph/9701387 . 



9. J. Kapusta, 'Finite- Temperature Field Theory', Eds. Landshoff et al., Camb. 1989, 
and references therein. 

10. J. Berges and K. Rajagopal, |hep-ph/9804233| ; M.A. Halasz et al., |hep-ph/9804290| . 



11. M. Stephanov, Phys. Rev. Lett. 76 (1996) 4472. 

12. R.A. Janik, M.A. Nowak, G. Papp and I. Zahed, Phys. Rev. Lett. 77 (1996) 4876. 



